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ABSTRACT 

Over 500 orbits of the synchronous satellites Syncom 2 , Syncom 3 
and "Early Bird" have been examined to reveal the long term effects of 
the longitude dependence of the earth's gravitational field. The period of 
record extends from August 1963 to January 1967. The east-west drift 
of these synchronous orbits are strongly influenced by the low order 
gravity harmonics and which are in resonance with them. The 
record shows these harmonics to have the following amplitudes and phase 

angles: = 1.80 ± .04, -15.3 ± 0.6°, 10®J = 0.16 

^33 = 26.4 _g* 3 Q. The record also shows, but less significantly, the in- 
fluence of the resonant harmonic H 3 j . The amplitude and phase of this 

+2 6 "i"73 ^ 

harmonic is found to be: 10 ®J 3 j = 1.5 ^ 3 ^ = -79 The bounds 

on these harmonic parameters are felt to be quite conservative and re- 
flect likely effects of higher order resonant gravity harmonics which 
could not be well determined by the data. The full record through Jan- 
uary 1967 now permits the discrimination of the long term east-west 
drift acceleration of a geostationary satellite to within 5.5 x lO"® deg/day^ 
(Maximum la) around the equator. 
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RESONANT GRAVITY HARMONICS FROM 3-1/2 YEARS OF TRACKING 
DATA ON THREE 24 HOUR SATELLITES 


INTRODUCTION 

Since July 1963 the National Aeronautics and Space Administration (NASA) 
and the Communications Satellite Corporation (COMSAT), and the Department 
of Defense (DOD) have tracked and operated eight synchronous (24 hour) satellites 
positioned at a fair number of longitudes around the equator. The free drift 
tracking record of three of these satellites (NASA's Syncoms 2 and 3 and COM- 
SAT's "Early Bird") prior to July 1965 has already revealed new and accurate 
information about the longitude dependence of the earth's gravity field^’^. The 
lowest or 2nd order longitude components of the field dominate the libratory 
east-west gravity drift of 24 hour satellites and can cause accelerations as 
high as 1.9 x 10'^ degrees per day.^ In order to be able to predict accurately 
their behavior at all longitudes, good tracking information must be obtained over 
as wide a longitude range as possible. This is true principally because third 
and fourth order components of the longitude gravity field have small but significant 
long term effects on these satellites. These effects can only be accounted for 
adequately by a wide longitude survey of synchronous satellite drift. This report 
brings the available tracking record of Syncom 3 and Early Bird up to January 
1967. It covers completely the longitude range of these satellites prior to that 
date. Complete tracking data on Syncom 2 from July 1963 up to April 1966 is also 
recorded here. Data on this satellite during the rest of 1966 adds about 6 de- 
grees of longitude to this survey. It will be incorporated in future reports on 24 
hour satellite drift which will also include the record of the NASA and COMSAT 
24 hour satellites ("INTELSAT" and "ATS") launched since December 1966. 

The tracking data presented in this report is analyzed for long term east- 
west drift acceleration according to the methods developed in previous investi- 
gations by the author^’'*. The measured accelerations are then used to solve 
for specific longitude dependent gravity terms in the earth's potential which can 
give rise to these perturbations. The essential point of this analysis is that the 
critical longitude gravity forces are in librational resonance on these satellites. 
The libration period due to the strongest effect is not less than 2 years. Thus 
the detailed tracking information taken over a few days, that is summarized by 
a set of 6 orbital parameters , is probably not much more valuable than the 6 
elements themselves in determining the disturbing field. Whenever small 
perturbation forces give rise to large long term effects, such as the orbital 
regression and precession due to the oblate earth, the most economical analysis 
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for those forces deals directly with the derived orbital elements themselves as 
the observed data®’®. The satellite's equator crossing longitude, while not a 
typical orbital element, is easily derived from the set of 6 that are generally given. 

The operators of 24 hour satellites are concerned mostly with the drift 
accelerations which can be expected around the equator. These are finally esti- 
mated through the disturbing earth potential derived from the set of observed drift 
accelerations on the synchronous satellites. The ultimate aim of this report as 
in the previous one”*, is to give realistic estimates of both the total acceleration 
and the individual gravity terms affecting the 24 hour satellites. 


ANALYSIS AND DATA REDUCTION 

The Earth's gravitational potential which is employed in the data analysis 
is the familiar spherical harmonic expansion in associated Legendre functions 
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(Sin T) 


Cosm/V + S^ij^ SinmX} , 


( 1 ) 


where /x is the earth's Gaussian gravity constant (3.98601 x 10® km®/sec^), r 
is the radius from the earth's center of mass to the satellite's, is the mean 
equatorial radius of the earth (6378.16 km), cp is the satellite's geocentric 
latitude and A. is its geocentric longitude The gravity coefficients and 
represent longitude dependent harmonic terms (H^J whenm^ 0. These are 
the terms which concern us here, and in particular, for the high altitude 24 hour 
satellites, only the lowest orders (-t) of these terms (approximately for i S A), 
because of the distance dumping factor (r^/r)"^ in the potential function of equa- 
tion (1). The more easily visualized amplitudes and phase angles of the 
harmonic terms (H^ J are related to the C^, coefficients of these terms by 
the formulas: 


=^tan-‘ 
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With these coefficients, the general disturbing harmonic term Ho in Equation 
(1) can be written as: 


= (7 ) fic. J^.Cosm (3) 

It is seen from Equation (3) that the longitude harmonic ^ through 

m imdulations in 360° of longitude. Thus we can speak of their wave length as 
360°/m and their longitude frequency as m (more accurately, m is the wave 
number: the number of longitude waves in over 360°). 

The long term east-west drift of the 24 hour satellite is dominated by a 
resonance with the term of the earth's potential which has the characteris- 
tics of pendulum motion^’® . In fact all the longitude dependent terms in Equation 
(1), contribute to this resonance as long as ^ - m is even. However, due to the 
distance damping of the anomalous potential, the term dominates the drift 
evolution. The dynamics of this long term pendulum-resonant drift is summarized 
by the following coupled equations^, giving the motion of the equator crossing 
longitudes (^ ) and the semimajor axis (a) of the orbit: 


and 
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The relevant coefficients Fo , to fourth order, ^ < 4, are given as: 


F 22 = [(Cos i + l)/2] 
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(6) 


The synchronous semimajor axis [in units of earth radii in Equation (6)] , is 
that semimajor axis at which k = 0. The orbital inclination is i . 

Equations (4) and (5) can be combined to eliminate their dependence on the 
forcing terms and yield an equivalent of Kepler’s period law. Thus: 


\ (- 4773^/1277 2) = a, or 

(7) 

A. (1/3t7) + a/a^ = 0. 


Integrating Equation (7) and, evaluating the constant of integration by specifying 
the synchronous condition, A. = 0 when a = a^ , yields: 


A. = - 3t 7 (a - a^)/a 


S 


rad. 

sid. day 


( 8 ) 


A first integral of Equation (4) can be immediately written down upon the separa- 
tion of the variables (A.)^ and \ , since = d (\)2/2dA., yielding: 


(/)2=c - 24772 (F^^/m) Cosm (\-/^^), (rad./sid. day)' 

■i-m EVEN 


( 9 ) 


where C is an arbitrary constant of integration. 
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Equations (4), (5), (8) and (9) constitute the principal condition equations 
governing the evolution in longitude and semimajor axis, of the 24 hour satellite 
due to the earth’s anomalous potential. Using position (\) and velocity or semi- 
major axis (a) data derived from the tracking of these satellites, we can solve 
these equations for the relevant longitude gravity harmonics which govern the 
drift. The semimajor axis data is taken directly from the sets of published mean 
elements for these satellites (Appendix A). The equator crossing longitudes (X) 
are derived from the published sets (Appendix A) of vector (or osculating) ele- 
ments, by numerical integration for less than one revolution in a realistic earth- 
moon-sun gravity field. 

Where the satellite’s drift rate is slow ( < .15°/day) and the crossing longi- 
tudes do not change significantly, the longitude acceleration and semimajor axis 
drift rate are essentially constant. For these drift arcs a low degree polynomial 
in the time can be used to determine the acceleration from the crossing or semi- 
major axis data directly. Where the satellite’s drift rate is more rapid it has 
been found to be most convenient to analyze the data according to the energy 
integral of the perturbed motion. Equation (9), Consecutive observations of the 
crossing longitudes from the same or neighboring sets of elements can serve to 
define the longitude rate over a small longitude range. Alternately, if semimajor 
axis data is given. Equation (8) defines the drift rate, with the mean synchronous 
semimajor axis being approximately 6.611 earth radii for both the 33° inclined 
and equatorial orbit satellites. 

For the purposes of this report, all the data in the fast drift arcs were re- 
duced to a single best determined longitude acceleration near the middle of these 
arcs. This acceleration was determined from Equation (4) using resonant gravity 
coefficients determined by a least-squares fit of semimajor axis or drift rate 
data according to Equations (8) and (9). In order to achieve the greatest precision 
in this acceleration measurement, with the limited data in each arc, only the 
dominant gravity terms (H 22 and occasionally were solved for in this pre- 
liminary fit. There is some gravity information lost in this preliminary smooth- 
ing process for the fast drift arcs. This will be recovered in future analyses of 
the slow drift (acceleration) data, and the fast drift (drift rate) data combined in 
a simultaneous solution for all the relevant gravity harmonics. 

In Appendix A is found the sets of elements and related tracking data in 13 
separate free drift 24 hour satellite arcs. Some of this data, and its reduction 
to acceleration-points, has already been reported”*, namely data arcs 1, 2, 3, 4, 

5, 6, 7, 8-1 and 9-1. A summary of this acceleration-point data reduction (proc- 
essed by the methods described above for the slow and fast drift arcs) is pre- 
sented for all the arcs, old and new, in Table 1. Only the old data arc 3 was re- 
jected as too poorly determined to be useful in the gravity solutions. The new 
drift arcs (10, 11, 12 and 13) of Syncom 3 and Early Bird came about from the 
following orbit maneuvers: 
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Arc 10, Syncom 3 (March- May 1965). 


The last orbit maneuvers on Syncom 3, performed by NASA, took place in 
the latter half of March 1965, resulting in a mean motion change from a slow 
westward to a slow eastward drift at about 172° east longitude. Subsequently 
this satellite was tracked and operated by DOD (Air Force System Command, 
Sunnyvale, California). From the last week in March to the end of April there 
was free drift from 172° to 173 1/2°. An orbit maneuver, only affecting the mean 
motion, reversed the slow eastward drift to a slow westward drift at the end of 
April. Free drift followed till the satellite reached about 172° at the end of 
May 1965. The total period from the end of March till the end of May has been 
analyzed as a single, slow, free drift arc (10) with a break in the drift rate on 
26 April 1965 (Appendix A). An orbit maneuver on Syncom 3, apparently took 
place in early June 1965 and again on 16 July. Free drift was allowed between 
16 July and 5 October 1965, but there are too few observations available during 
this period to compute a good acceleration. 


Arc 11, S 3 Ticom 3 (October — March 1965/66) 

Orbit maneuvers on Syncom 3 were performed on 5, 19 and 29 October 1965. 
From 29 October to June 1965/66 Syncom 3 was in a period of slow free drift 
from 172° to 165-1/2° east longitude. The complete record through March 1966 
is analyzed from the data in Appendix A. 


Arc 9-2 to 9-5 , Early Bird, (July - Nov. 1965) 

The Early Bird Arc labeled 9 (April - June 1965) in Reference [4] is here 
called 9-1. The free drift period for Early Bird between the end of April and the 
end of 1965 is actually only broken by a major mean motion change maneuver on 
2 December 1965. Four newly analyzed portions of free drift Arc 9 (Early Bird) 
are found in Appendix A, extending from 28° to 38° west longitude. 


Arc 12, Early Bird (December - May 1965/66) 

After the maneuver on 2 December 1965, the next Early Bird maneuver 
occurred on 20 September 1966. Five sub-arcs in this new free drift period for 
Early Bird have been analyzed in Appendix A. They cover the available record 
to May 1966, and extend from 38° to 28° west longitude. 
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Arc 13, Syncom 3 (September — January 1966/67) 


Orbit maneuvers were performed on Syncom 3 in June and late August 1966. 

A long, very slow, free drift period ensued from September 1966 through January 
1967, carrying the satellite from 161° to 160° east longitude. 

In addition to these new free drift arcs, the arc labeled 8 (Syncom 2) in 
Reference [ 4 ], has been extended to cover the continuing slow free drift of 
Syncom 2 from May 1965 to March 1966 between 65° and 83° east longitude. 

The estimates of model bias in Table 1 require some explanation. Initially, 
the crossing and semimajor axis data are reduced to acceleration points accord- 
ing to the slow or fast drift models discussed previously. Many numerical calcu- 
lations have shown^ that sun and moon gravity effects on 24 hour satellites cause 
long term east-west accelerations of about 0.03 x 10“® radians/sidereal 
day^ (Icr). Calculations also show that radiation pressure without consideration 
of the earth's shadow, causes a long term acceleration as high as 0.006 x 10'® 
rad./sid. day^ on the Syncom-Early Bird satellites. Taking into account the 
earth's shadow, I have estimated (from numerically calculated trajectories), the 
radiation pressure effect to be 0.086 x 10“ ® rad./sid. day^ (lo"). For many of 
the data arcs, more exact estimates of the gravity biases are available. These 
biases, due to the sun and moon as well as an incomplete earth model (in the 
fast drift reductions) were calculated from numerical trajectories closely 
paralleling the actual trajectories in these arcs. These calculated biases are 
likely to have errors of 0.005 x 10'® rad./days^ (la) due to short period 
effects. The numerical trajectories were calculated by Goddard Space Flight 
Center's "ITEM" Program (Interplanetary Trajectory by an Encke Method) in 
the presence of the sun and moon as well as an earth with realistic 2nd and 3rd 
order longitude gravity effects. For the fast drift arcs reduced to accelerations 
with only a 2nd order longitude gravity model (including the effects of only = 
Hjj for m 0), the total model biases are estimated to be between 0.040 and 
0.045 X 10"® rad./ sid. day^(lc7). For the slow drift arcs, the total model 
biases are estimated to be 0.035 x lo'® radians/sid. day^ (Icr). Where the 
long term gravity model biases are calculated, the remaining bias, from radia- 
tion pressure and short period gravity effects, is estimated to be 0.010 x 10"® 
( 1 ^). 

Measured accelerations, due just to the earth’s longitude gravity field, are 
calculated by adding these acceleration bias estimates to the acceleration re- 
duced observations. In Table 1, two sets of measured accelerations are tabulated: 
one with the calculated major sun-moon bias removed, and one without this ad- 
justment. In both sets, the bias deviations have been added to the measurement 
deviations. Total lo- measurement errors are estimated to be the square root 
of the sums of these deviations. 


7 


RESULTS 


Table 2 summarizes the gravity results of fitting the acceleration point data 
in Table 1 to the resonant drift model of equation (4). Of the five gravity terms 
whose maximum effects should be discernible by this data^, (Hjj, 
and ), only four had significant effect in reducing the 

residuals (s ) of the experiment. Least squares fits of both the adjusted and 
unadjusted (unweighted and weightec^ data to combinations of these coefficients 
show a strong reduction in the residuals (over a simple fit) only when the 
terms of H 33 are included in the model. A comparison of Test 1 with 2,3, and 4 
shows this strong reduction in both the residuals and the standard deviation for 
the dominant coefficients. The validity of the solutions in these tests, 
is also strengthened by the consistency of the coefficients shown in all the tests 
except perhaps in those for a full field to fourth order (excluding Even in 

the full field solutions (tests 9 and 10), the divergent values of compen- 

sated by a considerably larger standard error in the coefficients which is suffi- 
cient to overlap the simpler field solutions. With good data having sufficient 
longitude coverage, we would expect any of the terms H 3 J , H 33 and to 
show consistent solutions regardless of the model assumed. Even the small 
effect of H ^2 separable from which has the same longitude frequency, 

when better data on the moderately inclined Syncom 2 becomes available. As 
shown in Reference [4] , has almost no effect on Syncom 2, but a discernible 
one on the equatorial satellites. The difference would be detectible, resulting in 
a significant solution for > if enough good acceleration data on both Syncom 2 
and the equatorial satellites [with errors less than 0.015 x 10 rad./sid. day ^ 
(Icr)] were known. 

It is apparent from tests 1 - 10 in Table 2 that the strongest solutions for the 
harmonics, with this limited record, are achieved with the H 22 , H 33 combination 
tests (2,3, and 4) even though including 113 ^ and in the solutions (Tests 9 
and 10) result in the lowest residuals. The harmonic coefficients themselves 
(H 22 and H 33 ) are best determined in Test 3 on unadjusted but weighted data. 

The residuals in this test are about 0.013 x 10"® rad./sid. day^ greater than 
in the test with the least residuals (10). But the major coefficients appear better 
determined in the simpler field test because the solution overadjusts unrealistically 
for the minor effects in the complex field solutions. 

Comparing the tests on weighted unadjusted data (3, 5, 7 and 9) against those 
on weighted adjusted data (4, 6 , 8 and 10), we see apparently the same kind of 
unrealistic overadjustment for small effects damaging the solution with the more 
accurate data. The unadjusted data solutions are not heavily constrained towards 
critical data points. They are generally smoother solutions with somewhat 
smaller coefficient variance. However, the coefficient variance differences are 
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minor and the adjusted data solutions are more faithful to better data. Therefore 
I prefer these solutions in all cases where more than Hjj and H33 are being 
solved for and it is important to be as sensitive as possible to the best measure- 
ments. The tests for H22 and H33 (2,3, and 4) show little to choose between the 
use of adjusted, unadjusted, weighted or unweighted data. In both adjusted and 
unadjusted data tests, there is a noticeable reduction in the residuals when H3J 
is added to the model but at a small sacrifice in the H22 and H33 coefficient 
variances (compare tests 3 and 5, and 4 and 6). The smallest residuals of all 
are achieved with the adjusted data and a four term fit (Test 10). However, the 
variances in H22 appear to be excessive in this test. I believe the longitude 
survey is still too limited in quality and extent to provide realistic four term 
estimates. 

The lowest residual three term test with reasonable variances appears to be 
Test 6, and this was finally chosen as the basic solution. This solution is com- 
pared point for point with the adjusted measured data, in Table 3 (Column 1 vs. 
column 2). It is seen in this comparison that all points in the solution are within 
3o- of the measured data, 97% of the points are within 2 a- and 75% of the points 
are within la. The same number of points have positive as have negative resid- 
uals. These overall statistics are quite favorable. If the residuals were normally 
distributed, more than 99.7% of them should be within 3a, 95.5% within 2a and 
68.5% within la. The test 6 solution is a bit better than all the expectations 
which, of course, are based on an infinite number of data. 

The gravity parameters and covariance statistics of the basic solution is 
also illustrated in Figure 1, which is a calculation, from Equation (4), of the 
east-west acceleration on a geostationary satellite around the equator. The 
standard error estimation of this acceleration (the 10a curve) is particularly 
revealing. It conforms quite well with the measured data statistics in Table 3. 

A representative sample of these statistics are displayed as solid vertical lines 
in Figure 1. The overall 10a curve of the basic solution clearly reveals those 
longitudes where data is most needed to strengthen the resonant gravity solution 
for 24 hour satellites. The overall fit of the basic solution shows ala variation 
which is always less than 0.1 x 10 rad./ sid. day 2. Arbitrarily, I use an error 
of 0.075 X 10"® rad./day2 as defining the upper limit of acceptability for the 
solution. This defines three longitude zones where new data of reasonable quality 
can be expected to have the greatest effect in strengthening the overall validity 
of the gravity solution. Positioning a 24 hour satellite between about 12° and 55° 
east longitude, where no data exists, would probably yield the greatest improve- 
ment in geodetic information. Satellites in the zones between about 100° and 
157° west longitude, and between approximately 100° and 119° east longitude 
would also strengthen the 24 hour gravity solution considerably. (These con- 
clusions are also supported in Reference [10] .) By a strengthened solution, 
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1 mean most of all the reduction in the rather large uncertainty still remaining 
in the determination of the third order resonant gravity effects. Figure 2 shows 
the magnitude of this problem. From the error statistics of the basic solution, 
lo probability ellipses (solid curves) have been drawn in the H22, and H33 
quadrants enclosing areas of likely solutions with the limited model as seen by 
the data (See Appendix B). The ellipses in the H22 and quadrants are quite 
extensive, especially the ellipse, compared with that defining the variability 
of the H22 effect. Perhaps more important is the large amount of cross correla- 
tion in the basic solution between many of the gravity coefficients. This indicates 
most surely that effects of 2nd and 3rd order are not yet well separated. Once 
again, this is due undoubtedly to the limited extent of the longitudes surveyed 

so far by 24 hour satellites. 

As stated in the introduction the aim of this investigation is to derive 
realistic values for the individual gravity terms, as well as good estimates for 
the east-west accelerations. Therefore it is necessary to consider the likely 
variance in the basic solution due to the neglect of higher order resonant gravity 
effects which cannot be well solved for from the limited data available. In fact, 
solutions 2, 3, and 4 compared to the basic solution may be considered to be an 
example of such an estimate, since Hjj is obviously not well determined in the 
basic solution. As displayed in Figure 2 (x symbol) solution 4 gives estimates 
of H22 and H33 outside the basic solution Icr ellipses for these harmonics. 
Solutions 7, 8, 9 and 10 with other three and four parameter models also illustrate 
additional variance possibilities in H22. H3J and H33 . We can also test for 
limited-model bias by assuming unmodeled effects here have been well determined 
by other geodetic satellites- and ground surveys. Table 4 lists coefficient values 
from a number of such recent determinations. All of these except Izsak’s 1965 
values (which were updated in the SAO 1966 determination) are displayed in Figure 

2 for comparison with the 24 hour satellite results. Fixing 113^ , 

according to the results of Kaula (1966a) as an example. Solution 4 (Table 2) with 
the limited H22, H33 model (adjusted-weighted data) becomes Solution 12). 

Fixing H^2 by Kaula's 1966a values, the basic solution (6) becomes 

Solution 11. It is interesting, and perhaps significant that Solution 11 is little 
different from the basic solution in spite of the inclusion of higher order effects 
from an external source. On the other hand Solution 12 shows a rather strong 
divergance (and increased residuals) from Solution 4. It is seen from Figure 2 
that there is quite a large separation of the Solution 6 113^ determination from 
those of recent investigations. We know from the effect on the residuals that 
H3^ does have a noticeable influence on the solution. The discrepancy in 113^ 
as seen by the 24 hour satellites, while not as great as with still more limited 
data [compare with the symbols in Figure 2] is still sufficient to suggest there 
may be more uncertainty in both the and H33 results than is shown in the 
statistics of the basic solution. This conclusion is supported by the previous 
comparison between the simple two parameter solutions (2,3, and 4) and the 
basic three parameter solution (6). 
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Summarizing, the dashed rectangles in Figure 2 are realistic estimates of 
the likely absolute bounds on and H 33 as assessed by the 24 hour 

satellites so far (see Appendix B). They take into consideration the model simpli- 
fications forced on the solution due to the limited data available through 1966. 

The total variability in the coefficients from the basic solution over all the tests 
in Table 3, was the main criteria in setting these realistic bounds. Some dis- 
cretion was exercised, however, in crediting the extremes of the variability in 
Hj 2 from the simplest solution. Test 1, and in H 22 and Hjj from the most 
complex andprobably overadjusted solutions. Tests 9 and 10 . 

An interesting aspect of the results of the statistics of the "basic solution" 
is that the Icr ellipse does point in the direction of the most recent geodetic 
determinations (Figure 2). In fact, the 2 a- Hjj ellipse actually overtakes these 
values, but without enclosing them by a small margin. However, an examination 
of the correlation coefficients of this solution is not as encouraging as this 
simple observation suggests. It is true that the coefficients are strongly 

and positively correlated. Thus movement in the ellipse towards the pre- 
dominant recent determinations is also correlated to movement in the H 33 ellipse 
towards these same solutions. However, the Hgj correlation coefficients 
and C 33 are strongly negative so that the movement towards the predominant 
recent solutions in the and H 33 ellipses is strongly correlated to movement 
away from these solutions in the H 3 j ellipse. Evidently, the direction of these 
ellipses suggest better agreement to recent results than is actually indicated 
by the data. 

A further comparison of three of these recent determinations, directly with 
the adjusted 24 hour satellite measurements, is displayed in Table 3. From this 
comparison, it is apparent that the solution which best represents the 24 hour 
observations is that of Kaula (1966a). The Kaula 1966b solution (A symbols in 
Figure 2 ) would, evidently, compare best of all with the 24 hour data, but that is 
because this solution incorporates in it a fair amount of this data. In fact the 
Kaula (1966a) solution, indirectly is also weighted towards the 24 hour data; but 
to a smaller degree than that of Kaula (1966b). Considering independent solutions 
along, SAG (1966) appears to be a somewhat superior solution than Anderle (1965) 
on the basis of the 24 hour satellite data. Comparing the two Doppler-Satellite 
solutions from Figure 2 , the Anderle one would probably prove somewhat better 
overall, than that of Guier and Newton, because of its closer H 22 agreement with 
the 24 hour satellite solution. But clearly, all of these recent solutions have 
failed to predict even reasonably well, the accurate accelerations measured in 
the Early Bird arcs. 

A graphical display of the unadjusted and adjusted and predicted accelerations 
in the Early Bird and Syncom 3 arcs as well as in the last Syncom 2 drift period, 
is shown in Figure 3. 
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As incidental, but important statistics from the basic solution, the four 
equilibrium points (zero acceleration) for the geostationary satellite are located 
at (See Figure 1): 

X.J = 75.9 ± 0.6° (Stable Equilibrium) 

^2 = 161.9 ± 0.3° (Unstable Equilibrium) 

= -106.7 ± 0.9° (Stable Equilibrium) 

= -11.9 ±0.5° (Unstable Equilibrium). 

All of these values except for are within the lo- estimates given in Reference 
[4] from more limited data. The bounds on k^ given here overlap those de- 
termined in the earlier study. 

The maximum east-west acceleration on the geostationary satellite, due to 
the earth, occurs at 119° east with a magnitude of -3.18 ± 0.08 x 10“® rad./sid. 
day^ = -1.83 ± 0.5 x 10’^ deg./solar day^. This is identical to the result of 
the earlier study'*, so that, to correct continuously for this east-west acceleration 
would require (as calculated previously) conservatively, a velocity increment of: 


AV = 6.38 ft./sec./yr. 


Finally, as a check on these results, we can compare the coefficients seen 
in this new solution, with those seen in the earlier study'*. The new solution 
(Table 2 and Figure 2) is, except for the ill determined coefficients, within 
the Icr bounds found in the earlier study. The coefficients absolute bounds in 
this new solution are, however, greater than those in the old study, in spite of 
the fact that more and better data has been used here. This conservatism is 
due, principally, to the method of evaluating the model bias effects. Here these 
effects are estimated by direct comparison of many parameter tests on the actual 
acceleration record, with and without use of external gravity information. In the 
older study, the likely model bias was evaluated only indirectly from simulated 
data which was of a much higher quality than the actual acceleration record. 

As a result, it appears from this study that only H 22 (rather than both H 22 and 
H 32 ) is significantly better determined from the available 24 hour satellite record, 
than in any of the other recent geodetic reductions. 


12 


/ 


CONCLUSIONS 

1. As in the previous study of 2 years of the 24 hour satellite record^, I 
conclude that virtually all of the east-west geographic acceleration of 
the 24 hour satellites can be accounted for by the second and third order 
sectorial harmonics of the earth's gravitational field which resonate with 
them. 

2. With due adjustment for the small long term effects of the sun and moon's 
gravity and radiation pressure, and considering the neglect of likely higher 
order resonant earth gravity, these dominant sectorial harmonics are 
estimated to be: 


10® J 22 = 1.80 ± 0.04, 

which corresponds to a difference in major and minor radii of 68.8 ± 1.8 
meters in the earth's elliptical equator, and 


\22 = -15.3 ± 0.6°, 

which is the longitude location of the major axis of the elliptical equator, 
and 


J33 


0 

^•^°- 0.03 




33 


26.4+®-®° . 
-8.3° 


3. The sectorial harmonics above, within their ranges, are believed to be 
absolute or true measures of those individual components of the. earth's 
field. 

4. A third resonant earth harmonic, , was evident but poorly discrimi- 
nated from the limited acceleration record. The data shows the absolute 
bounds on this harmonic to be: 
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10® J31 = , and 




31 


= -79 


+73° 
- 88 ° ■ 


5. The 3-1/2 year record of 24 hour satellite accelerations shows that the 
geostationary satellite can be in uncontrolled long term east-west 
equilibrium at only the following four longitude locations; 

= 75.9±.6° (Stable Equilibrium) 

= 161.9 ± .3° (Unstable Equilibrium) 

= -106.7 ± .9° (Stable Equilibrium) 

= -11.9 ±. 5° (Unstable Equilibrium). 

6. The maximum long term longitude acceleration on the geostationary 
satellite is, as reported previously: 

^ = -1.88x 10“^ deg./day2, 

occurring at about 119° east of Greenwich, and requiring a velocity 
increment of 6.38 ft/sec/year to counteract continuously. 

7. The H22 term, measured from this 3-1/2 year data record, is probably 

the most accurate known at this time. However, the accurate measure- 
ment of H33 and higher order effects must wait on receipt of new 

high quality tracking data, especially from longitudes not covered by the 
operating 24 hour satellites studied here. 
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Figure 3A. Drift Acceleration in Syncom 2 Arc 8. 
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ncom 3 Arcs 6, 7, 10, 11 and 13. 
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Figure 3C. Drift Acceleration in Early Bird Arcs 9 and 12. 



Table 1 

Data From 13 Free Drift Arcs of Three 24 Hour Satellites 


(D® 


SATELLITE 

ARC 

LONGITUDE 

{^): 

DEGREES 

— 

SEMIMAJCNI 
AXIS (a): 
EARTH RADD 

INCLINATION 

uy. 

DEGREES 

MEASURED DRIFT 
ACCELERATION 
(V)-. 10 ® RADIANS 
PER SIDEREAL 
DAY^ 

" (V), 

MEASURED: 

10 * RADIANS 
PER SIDEREAL 
DAY 2 

« ( ) BIAS, 
ESTIMATED: 

10 * radAiay^ 

^ BIAS, 
CALCULATED: 
10 * RAD/DAY2 

X. MEASURED 
ACCELERATION, 
UNADJUSTED FOR 
CALCULATED BIAS: 
10 -* R/D^ 

X. MEASURED 
ACCELERATION, 
ADJUSTED FOR 
CALCULATED 
BIAS: 10 * R/D^ 

SYNCCM 2 

1 

-55.22 

6.6111 

33.02 

-2. 253 

0.033 

0.035 

± 0. 015 

±0.010 

-2.253 

±0.048 

-2,238 

±0.035 

SYNCOM 2 

U 

-55.235 

6. 611 

33.02 

2.255 

0.078 

0.035 



-2.255 

±0. 086 

-2.255 

±0.086 

SYNCOM 2 

2 

-60.94 

6.6110 

32.825 

-2.291 

0.057 

0.035 

-0.029 

±0.010 

-2.291 

±0.067 

-2.320 

±0.058 

SYNCCaa 2 

2a 

-61.00 

6.012 

32.83 

-2.268 

0. 117 

0.035 



-2.288 

±0.122 

-2.288 

+0 122 

SYNCOM 2 

4 

-140.00 

6.6204 

32. 58 

2. 136 

0.084 

0.040 

+ 0. 015 

±0.010 

2. 138 

±0.093 

2.153 

±0.t)65 

SYNCCM 2 

4a 

-140.00 

6.620 

32.58 

2.386 

0. 155 

0.040 



2.366 

±0. 160 

2.366 

±0.100 

SYNCOM 2 

5-1 

161.00 

6.6165 

32.40 

-0. 199 

0.066 

0.040 

+ 0.005 

±0.010 

- 0. 199 

±0.077 

- 0. 194 

±0.067 

SYNCOM 2 

5a 

130.00 

6.017 

32.3 

-2. 550 

0.082 

0.045 



- 2. 550 

±0.093 

- 2. 550 

±0.093 

SYNCOM 2 

5-2 

104.50 

0.6178 

32.15 

-2.278 

0.068 

0.040 

- 0. 041 

±0.010 

-2,278 

±0,077 

-2.319 

±0.067 

SYNCOM 3 

e 

178.71 

6. 6115 

0.11 

1.707 

0.059 

0.035 

-0.024 

±0.010 

1.707 

±0. 069 

1.683 

±0.060 

SYNCOM 3 

7 

176.80 

6. 6123 

0. 27 

1.550 

0. 175 

0.035 

+ 0.027 

±0.010 

1.550 

±0. 178 

1.577 

±0. 175 

SYNCOM 2 

8-1 

06. 115 

6.6112 

31.89 

0.950 

0. 062 

0.035 

+ 0.018 

±0.010 

0.950 

±0.071 

0.968 

±0.063 

SYNCOM 2 

8- la 

66.39 

6.61125 

31.89 

0.845 

0.078 

0.035 



0,845 

±0.086 

0.845 

±0. 086 

SYNCOM 2 

8-2 

65.31 

6. 81095 

31.79 

1.017 

0.040 

0.035 



1.017 

±0. 053 

1.017 

±0.053 

SYNCOM 2 

8- 2a 

65.30 

6.6109 

31.78 

1.000 

0. 040 

0.035 



1.000 

±0.053 

1.000 

±0.053 

SYNCOM 2 

8-3 

69.00 

0.61 

31.4 

0.676 

0.032 

0. 035 



0. 676 

±0.047 

0.670 

±0.047 

SYNCOM 2 

8-4 

77.80 

6.61 

31.4 

-0. 152 

0.021 

0.035 



-0.152 

±0.041 

-0. 152 

±0.041 

SYNCOM 2 

8-5 

73. 25 

6.61 

31.4 

0.235 

0.015 

0.035 



0.235 

±0.038 

0.235 

±0.038 

EARLY BIRD 

9-1 

-28. 70 

6.0105 

0.20 

-1.441 

0. 010 

0.035 

-0.031 

±0.010 

-1.441 

±0.036 

-1.472 

±0.014 

EARLY BIRD 

9-2 

-28. 37 

6.6112 

0.43 

-1.478 

Q.OlO 

0.035 

+ 0.030 

±0.010 

- 1. 478 

±0.036 

-1.448 

±0.014 

EARLY BIRDl 

9-3 

-28. 54 

6.0112 

0.43 

-1.453 

0.016 

0. 035 



-1.453 

±0.038 

-1.453 

±0. 036 

EARLY BIRD 

9-4 

-31.86 

6.61185 

0.55 

-1.682 

0. 014 

0.035 



- 1.682 

±0.038 

-1.682 

±0.038 

EARLY BIRD 

9-5 

-34. 57 

6.0123 

0.64 

-1.849 

0.008 

0.035 



-1.849 

±0. 036 

- 1. 849 

±0.038 

SYNCOM 3 

10 

172.75 

6.611 

0.00 

1,072 

0.093 

0.035 



1.072 

±0.099 

1.072 

±0.099 

SYNCOM 3 

10a* 

172.75 

6.611 

0.00 

0.907 

0.345 

0.035 



0.907 

±0. 345 

0.907 

±0. 345 

SYNCOM 3 

11-1 

169. 10 

6.61 

0.53 

0.831 

0. 018 

0.035 



0.831 

±0.039 

0.831 

±0. 039 

SYNCOM 3 

11-2 

168. 29 

6.61 

0.53 

0.699 

0.020 

0.035 



0.699 

±0.040 

^0.899 

±0.040 

EARLY BIRD 

12-1 

-35. 98 

6.6094 

0.74 

-1.954 

0.021 

0.035 

-0. 022 

±0.010 

-1.954 

±0.041 

-1.978 

±0.023 

EARLY BIRD 

12-2 

-31.46 

6.6101 

0.85 

- 1. 678 

0.010 

0.035 

+ 0.045 

±0.010 

- 1. 678 

±0.036 

-1.633 

±0.014 

EARLY BIRD 

12-3 

-29. 83 

6.6101 

0.85 

-1.556 

0,011 

0. 035 

+ 0,029 

±0.010 

- 1. 556 

±0, 037 

-1.527 

±0.015 

EARLY BIRD| 

12-4 

-31.29 

6. 6103 

0.90 

-1.686 

0.007 

0. 035 

+ 0.044 

±0.010 

-1.688 

±0.036 

- 1. 644 

±0.012 

EARLY BIRD 

12-5 

-28.83 

6. 6103 

0.90 

-1.467 

0.006 

0.035 

-0.027 

±0.010 

-1.467 

±0. 036 

- 1. 494 

±0.012 

SYNCCM 3 

13 

160.74 

6.61 

1.3 

-0,090 

0. 010 

0. 035 



-0.090 

±0.036 

-0.090 

±0.036 


NOTES <T)Da*A arcs 1, 2, 4, 5-1, 5-2. 6, 7. 8-1, AND »-l ARE IDENTICAL TO ARCS 1. 2. 4. 5A, 5-18, 6, 7, 8 AND 9 OF REFERENCE 1 4], TABLE 10. ARCS la. 2a, 4a. ETC. ARE 
ANALYZED FROM SEMIMAJOR AXIS DATA. SEE APPENDIX A FOR THE DATA USED IN ALL THE NEW ARCS. 

0 THESE ARE THE LONGITUDES NEAR THE MID POINTS OF THE DRIFT ARCS, WHERE THE ACCELERATION IS BEST DETERMINED BY THE DATA (SEE REFERENCE [4] 
AND APPENDDI A FOR THE LONGITUDE SPAN COVERED IN EACH DATA ARC. AND TABLE 3 FOR THE TIME SPAN IN EACH ARC). 

0 THIS IS THE ACCELERATION AS REDUCED FROM THE CROBSINC OR SEMIMAJOR DATA BY A SIMPIE LOW ORDER TIME POLYNOMIAL OR ENERGY ( DRIFT RATE ) 
INTEGRAL MODEL. AS DESCRIBED IN THE TEXT. 

0 THIS IS THE ESTIMATION OF THE MINIMUM STANDARD DEVIATION OF THE ACCELERATION, MEASURED AS DESCRIBED ABOVE AND IN THE TEXT. 

0THIS BIAS DEVIATION ESTIMATE IS BASED ON AN ESTOAATE OF *0.032 x lO'* RAD/ DA Y^ FOR SUN - MOON - EARTH MODEL BIAS AND ±0.008 x lO'* R/D^ FOR 
RADIATION PRESSURE, FOR THE SLOW DRIFT ARCS, AND AN INCREASE IN THE SUN- MOON- EARTH MODEL BIAS TO ±(0,040—0.045) x 10 -5 RAD/ DAY ^ FOR THE 
FAST DRIFT ARCS. 

0THIS IS THE SUN- MOON- EARTH MODEL BIAS CALCULATED FROM NUMERICALLY GENERATED TRAJECTORIES CLOSELY PARALLELING THE ACTUAL 
TRAJECTORIES (SEE REFERENCE (4J). THE UNCERTAINTY OF O.Ol x 10** RAD/DAY^ IS DUE TO UNCERTAINTY IN THIS CALCULATION PLUS LDCELV 
RADIATION PRESSURE BIAS. 

0 COLUMN 0 DATA WITH DEVIATIONS TAKEN AS THE SQUARE ROOT OF THE SUMS OF THE SQUARES OF THE DEVIATIONS IN COLUMNS 0 AND (s). 

0 COLUMNS® + ® DATA WITH DEVIATIONS TAKEN AS THE SQUARE BCWT OF THE SUMS OF THE SQUARES OF THE DEVIATIONS IN COLUMNS 0 
AND 0. 

♦ OMITTED FROM THE GRAVITY ANALYSIS; ACCELERATION TOO POORLY DETERMINED. 
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Table 2 

Solutions (Tests) for Gravity Coefficients From 24 Hour Satellite Data 
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Table 3 


Measured and Predicted Accelerations on 24 Hour Satellites 


® ® (D © © 


ARC 

SATELLITE 

LONGITUDE* 

(DEGREES) 

DRIFT TIME 

MEASURED** 
DRIFT ACCELERATION 
( 10-5 RADIANS /Sm. 
DAY^) 

PREDICTED DRIFT ACCELERATION ( 10’5 RAD/SID. DAY^) 

WAGNER (1967 b) 
(COEFnCIENTS OF 
TEST 6, TABLE 2) 

KAULA (1966 a) 

SAO (1966) 

ANDERLE (1965) 

1 

SYNCOM 2 

-55.22 

AUG -DEC *68 

-2.238 

A 0.035 

-2.234 

-2.146 2* 

-2.079^ 

- 2. 185 

la 

SYNCOM 2 

-55.24 

AUG -DEC ’63 

-2.255 

±0. 086 

-2.234 

-2. 146 

- 2 . 0792 * 

-2.185 

2 

SYNCOM 2 

-60.94 

DEC - MAR ’63 - ’64 

- 2. 320 

±0.058 

-2.254 

- 2 . 1392 * 

-2.069^ 

-2.200 2® 

2a 

SYNCOM 2 

-61.00 

DEC - MAR ’63 - ’64 

-2.288 

±0. 122 

-2.254 

- 2. 138 

-2.068 

- 2. 199 

4 

SYNCOM 2 

- 140. 00 

APR -JULY ’64 

2.153 

±0.085 

2. 189 

2.114 

2.066 

2. 140 

4a 

SYNCOM 2 

- 140. 00 

APR - JULY '64 

2.366 

±0. 159 

2.189 

2.114 

2.066 

2. 140 

5-1 

SYNCOM 2 

161.00 

JULY -FEB '64 -’65 

- 0. 194 

±0.067 

- 0. 141 

-0. 133 

-0.110 

-0. 004^ 

5a 

SYNCOM 2 

130.00 

JULY - FEB ’64 - ’65 

- 2. 550 

±0. 093 

- 2. 484 

-2.541 

-2.557 

-2.640 

5-2 

SYNCOM 2 

104. 50 

JULY -FEB ’64- *65 

-2.319 

±0. 067 

-2.326 

- 2. 433 

-2.438 

-2.623^ 

6 

SYNCOM 3 

178.71 

OCT- JAN ’64- ’65 

1. 683 

±0.060 

1.749 

1,668 

1.705 

1. 8812* 

7 

SYNCOM 3 

176.80 

JAN -MAR ’85 

1.577 

±0. 177 

1.575 

1.499 

1.536 

1.711 

8-1 

SYNCOM 2 

66. 12 

FEB -JUNE ’65 

0. 968 

±0. 063 

0.924 

0.993 

1.021 

0. 968 

8- la 

SYNCOM 2 

66. 39 

FEB -JUNE *65 

0. 845 

±0. 086 

0.900 

0. 967 

0. 995 

0. 940 

8-2 

SYNCOM 2 

65.31 

FEB -JUNE '85 

1.017 

±0.5s3 

0. 996 

1.069 

1.096 

1.048 

8-2a 

SYNCOM 2 

65.30 

FEB -JUNE ’65 

1.000 

±0. 053 

0.997 

1.070 

1.097 

1.049 

10 

SYNCOM 3 

172. 75 

MAR -MAY ’85 

1.072 

±0. 099 

1.161 

1.114 

1. 151 

1 . 3192 * 

8-3 

SYNCOM 2 

69.00 

JULY - MAR ’65 - ’66 

0. 676 

±0.047 

0. 666 

0.718 

0. 747 

0. 676 

8-4 

SYNCOM 2 

77.60 

JULY - MAR ’85 - ’66 

- 0. 152 

±0.041 

- 0. 160 

-0.161 

-0. 133 

-0.2522* 

8-5 

SYNCOM 2 

73. 25 

JULY - MAR *65 - ’66 

0.235 

±0. 038 

0.263 

0.289 

0.3182* 

0.223 

11-1 

SYNCOM 3 

169. 10 

NOV -MAY ’65 -’86 

0. 831 

±0. 039 

0.801 

0.740 2* 

0. 775 

0 . 9332 * 

11-2 

SYNCOM 3 

168.29 

NOV - MAY ’85 - ’66 

0. 699 

±0.040 

0.713 

0. 654 

0. 688 

0. 843 2* 

13 

SYNCOM 3 

160. 74 

SEPT - JAN ’86 - ’67 

-0. 090 

±0. 036 

- 0. 127 

-0.1812* 

-0.157 

-0.034 

9-1 

EARLY BIRD 

-28. 70 

APR -JUNE '65 

-1.472 

±0.014 

-1.463 

- 1.446 

- 1 . 4242 * 

-1.3612* 

9-2 

EARLY BIRD 

-28. 37 

JUNE - AUG ’65 

- 1. 448 

±0.014 

-1.437 

-1.423 

-1.4022* 

- 1 . 3372 * 

9-3 

EARLY BIRD 

-28. 54 

JUNE - AUG ’65 

- 1. 453 

±0.038 

-1.450 

-1.435 

- 1. 413 

- 1 . 3492 * 

9-4 

EARLY BIRD 

-31.86 

SEPT -OCT ’65 

- 1. 682 

±0.038 

-1.890 

- 1. 652 

-1.618 

- 1 . 5772 * 

9-5 

EARLY BIRD 

- 34. 57 

OCT -NOV '65 

- 1. 849 

±0. 036 

-1.869 

-1.811 

- 1.7692* 

- 1. 7462* 

12-1 

early BIRD 

-35.98 

DEC - JAN '65 - ’66 

-1.976 

±0.023 

- 1. 958 

- 1 . 8892 * 

-1.8422* 

-1.829^ 

12-2 

EARLY BIRD 

-31.46 

JAN - MAR ’66 

- 1, 833 

±0.014 

-1.6632* 

- 1.628 

- 1. 5962* 

-1.551»* 

12-3 

EARLY BIRD 

- 29. 83 

MAR - APR *66 

- 1. 527 

±0.015 

-1.546 

- 1. 522 

-1.4962* 

- 1. 4412* 

12-4 

EARLY BIRD 

-31. 29 

JAN - MAR '66 

- 1.644 

±0.012 

-1.651 

-I. 6 I 72 * 

-1.585^ 

- 1. 540*^* 

12-5 

early bird 

-28.83 

APR -MAY ’66 

- 1. 494 

±0.012 

-1,472 

- 1 . 4552 * 

- 1. 4322* 

-1.370®* 


* THESE ARE THE LONGTTTOES NEAR THE MID POINTS OF THE DRIFT ARCS, WHERE THE ACCELERATION IS BEST DETERMINED BY THE DATA (SEE 
reference [4] AND APPENDIX A FOR THE LONGITUDE SPAN COVERED IN EACH DATA ABC). 

’• ADJUSTED DATA FOR LIKELY SUN- MOON- EARTH MODEL BIAS (SEE TABLE 1). 
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Table 4 

Low Order Longitude Gravity Harmonics From Recent Investigations 
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APPENDIX A 
Basic Data 


In this appendix I list the basic data used in this gravity analysis. Except 
for Early Bird, the data consists mainly of sets of cartesian and spherical co- 
ordinate ("ADVARB") vectors, or mean Keplerian elements for the Syncom satel- 
lites 2 and 3, reduced from primary range and range rate and minitrack observa- 
tions. The Data Systems and Tracking Directorate at Goddard Space Flight Center 
(through the offices of Dr. Joseph Siry) supplied both cartesian and mean element 
sets and related equator crossing data (at near epoch times) for Syncoms 2 and 3 
from August 1963 through the spring of 1965. Subsequently, the Airforce Systems 
Command at Sunnyvale, California has been supplying "ADVARB" vector elements 
and near epoch time crossing data for Syncoms 2 and 3. 

The Early Bird data listed here consists mainly of equator crossing longi- 
tudes for the period April 1965 to May 1966, estimated directly from range, 
azimuth and elevation measurements near the crossings, taken from Comsat's 
tracking facility at Andover, Maine. 

Table Al, repeated for completeness from Reference [4], gives the GSFC 
vector elements for data Arcs 1 through 5. The analysis of these elements, re- 
sulting in the reduced accelerations listed for Arcs 1, 2, 4, 5-1, and 5-2 in Table 
1, has been adequately described in Reference [ 4 ]. (Arc 5-1 here was called 
Arc 5A in Ref. [4], and Arc 5-2 was called Arc 5-18.) Table A2 gives the 
same kind of data, already analyzed, for Arcs 6, 7 and 8. The results of the 
analysis of data Arc 8 in Table A2, is listed in Table 1 as Arc 8-1. 

In Tables A3 and A4 I list the mean semimajor axis and related data for 
Syncoms 2 and 3 , also reported by the GSFC Tracking Directorate which I have 
used to calculate additional, independently observed accelerations for data 
arcs: 8-la (slow drift) , 8-2 (slow drift) 8-2a (slow drift) , la (slow drift) , 2a 
(slow drift), 4a (fast drift), 6a (slow drift), and 7a (slow drift). The accelerations 
in Arcs 6a and 7a were too poorly determined to be useful in the gravity analysis 
and are not reported in Table 1. The best acceleration results for arcs 8-la, 

8-2, 8-2a, la, 2a and 4a are reported in Table 1. 

In Table A5 I list the mean semimajor axis and related GSFC reported data 
in the fast drift Arc 5 (Syncom 2), The single, best observed acceleration found 
from this data is reported as data Arc 5a in Table 1. 
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In Table A6 is listed tracking data for the "broken" arc 10 (Syncom 3). The 
best acceleration results for this slow drift arc are reported in Table 1 as data 
arcs 10 and 10a, though only the equator crossing analysis gave sufficiently small 
residuals to be useful for the gravity solutions. 

In Table A7 is listed a set of Department of Defense (DOD) elements, com- 
puted at Sunnyvale, California, from the spring of 1965 to January 1967. These 
vector elements (ADVARB) were computed from range and range rate tracking 
at Stations in Hawaii, the Phillipines, Vietnam and Africa. The earth model used 
to reduce the tracking data to these vector elements was more complex than that 
used to derive the GSFC elements. However, the accuracy of these elements 
when derived from information over a few revolutions, are not superior to the 
GSFC elements because of this, due to the very small short term effects of the 
more sophisticated model on the high altitude 24 hour satellites. Since the data 
used to derive the bulk of elements spans no more than a one week period, even 
the long term (resonance) effects of the longitude gravity model used, are rela- 
tively small. However, the earth model used by DOD, Sunnyvale, since the spring 
of 1966 1^ incorporated the gravity harmonics reported in Reference [ 4 ] and 
has resulted, on average, in better elements when the data span used is greater 
than about a week. 

These DOD elements have the following definition: Epoch, Date (Universal 
Time): M = month, D = day, Y = year, H = hour, MI = minute, S = second: a° = 
right ascension of the satellite with respect to the vernal equinox at epoch (de- 
grees), 5° = declination with respect to the equator at epoch (degrees), /3° = 
flight path angle from the satellite's outward geocentric radius vector to the 
velocity vector (degrees), AZ° = azimuth or the angle (from north clockwise) 
between the planes formed by the satellite's geocentric radius and the earth's 
north pole vectors and the satellite's radius and velocity vectors (degrees), RAD' 
= the satellite's distance (radius) to the center of mass of the earth (feet), V/FPS 
= the satellite's velocity in feet per second. The sets of DOD elements applicable 
to the long free arcs 8-3 to 8-5, 11, and 13 are shown in Tables A7 and A8. 
Additional information from these and other orbits, used in the long term accele- 
ration analysis for these slow drift arcs, are shown in Tables A9, AlO and All. 
The best acceleration results for these arcs are listed in Table 1. 

Finally, in Table A12 is listed the ascending and descending equator crossings 
for Early Bird in Arcs 9 and 12, as estimated directly from near-crossing range 
azimuth and elevation observations (supplied by the COMSAT Corporation) taken 
at the Andover, Maine tracking station. Also listed are occasional osculating 
(vector) Keplerian elements for Early Bird, determined from the full set of 
Andover Tracking data over a number of revolutions past the orbit date. The 
best acceleration results from these slow drift arcs are shown in Table 1. The 
data-arc time spans (as for all the acceleration points) appear in Table 3. 
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Table A1 


Inertial Position and Velocity Coordinates for Syncom 2 and 
Syncom 3 as Reported by GSFC* 


Tracking Epoch 

X 

Y 

Z 

X 

Y 

i 

(Yi-mo-day-ht-miB UT) 

(10< km) 

(10< km) 

(10* km) 

(km /sec) 

(km /sec) 

(km/sec) 


/63-8-18-1-30.0 

1.8517253 

-3.66 56 408 

-0.95346425 

2.5197630 

0.87570544 

1.5296730 


63-8-22-6-12.14 

3.8192813 

0.19653916 

1.7732339 

-0.64139671 

2.8111679 

1.0698498 


63-8-26-17.0 

-3.9190365 

0.71434463 

-1.3838910 

-0.030659063 

-2.7659652 

-1.3414586 


63-8-31.0 

1,2517473 

-3.8173186 

-1.2805028 

2.7082374 

0.42130417 

1.3937010 


63-9-3-13-23.0 

-2.6683093 

3.2430659 

0.37724660 

-2.0920997 

-1.5282086 

-1.6556391 


63-9-5.0 

1.5690433 

-3.7540418 

-1.1058861 

2.6179496 

0.66119258 

1.4710835 


63-9-9.0 

1.8101775 

-3,6842867 

-0,96169979 

2.5333080 

0.84735057 

1.5232783 


63-9-12-2.0 

3.4100141 

-2.4566117 

0.33394973 

1.4036497 

2.1745471 

1.6606540 


63-9-17-2.0 

3.5605032 

-2.1949356 

0.52762629 

1.1881299 

2.3199229 

1 .6320863 

Arc 1( 

63-9-20-2.0 

3.6381029 

-2.0310580 

0.64 240062 

1.0551433 

2.3993342 

1.6084890 


63-9-27-2.0 

3,7828113 

-1.6282253 

0.90343513 

0.73324914 

2.5577582 

1.5413518 


63-10-1-2.0 

3.8398851 

-1.3915057 

1.0449696 

0.54792106 

2.6322531 

1.4928190 


63-10-8-2,0 

3.8997221 

-0.96844341 

1.2765410 

0.22281006 

2,7325580 

1.3929604 


63-10-14-2.0 

3.9114166 

-0.59213158 

1.4569605 

-0.059266662 

2.7895298 

1.2936179 


63-10-22-2.0 

3.8655564 

-0.089237507 

1.6821569 

-0.43011588 

2.8231687 

1,1399466 


63-10-30.0 

3.7935711 

-1.5763511 

0.94831895 

0.69205550 

2.5804974 

1.5229227 


63-11-6.0 

3.8724273 

-1.1762815 

1.1818209 

0.38097748 

2.6938419 

1.4335987 


63-11-12-5.0 

1.1446010 

3.4554430 

2.1294725 

-2.7244144 

1.2840822 

-0.61733366 


k^3-ll-18-13.0 

-3.7038702 

-0.58064695 

-1.9332132 

0.90759960 

-2.7961257 

-0.89794623 


[43-11-28-1.0 

3.4781437 

1.1376687 

2.0944723 

-1.2933144 

2.7059520 

0.67780023 


63-12-4.0 

3.7151749 

0.53514140 

1.9201218 

-0.87415171 

2.8043813 

0.90960269 


83-12-10.0 

3.5731662 

0.92615699 

2.0379487 

-1.1461780 

2.7500585 

0.76033287 


63-12-16-17.0 

1.0227024 

-3.8691081 

-1.3292554 

2.7458357 

0.25931246 

1.3585847 


63-12-23-19.0 

3.0804533 

-2.8787714 

0,055186822 

1.7485346 

1.9038213 

1.6650503 


64-1-6-17.0 

2.2436953 

-3.5177574 

-0.60704656 

2.32626 1 5 

1.2062956 

1.6089899 


64-1-9-6.0 

-3.0926247 

2.8665793 

-0.080953346 

-1.7324029 

-1.9159081 

-1,6667974 

Arc a 

64-1-15-18.0 

3.3299438 

-2.5658898 

0.32807382 

1.4797710 

2.1321920 

1,6489369 


64-1-20-21.0 

3.7030072 

0.56000451 

1,9369923 

-0.88759725 

2.8080493 

0.88485319 


64-1-29-20.0 

3.8204264 

0.093177361 

1.7811029 

-0.55721471 

2.8372181 

1.0470801 


64-2-5-16.0 

2.8698222 

-3.0876991 

-0.10282631 

1 .9209748 

1.7308532 

1.6640514 


64-2-10-19.0 

3.8650196 

-0.19769561 

1.6738813 

-0.34675337 

2.8365031 

1.1356011 


64-2-17-17.0 

3.7268026 

-1.7536007 

0.90292256 

0.82362190 

2,5381818 

1.5281483 


64-2-25-19.0 

3.6533259 

0.69908427 

1.9851452 

-0,98120516 

2.7965023 

0.82109998 


64-3-4-23.0 

0.18113658 

3.8087583 

1.8019013 

-2.8347353 

0.61701878 

-1.0178939 


1^64-3-10-13.0 

2.0644779 

-3.6124046 

-0.68632246 

2.4052556 

1,0737921 

1.5862758 


[44-3-18-3.0 

-3.6249117 

2.0233288 

-0.74902788 

-1.0365741 

-2.4353607 

-1.5693664 


64-3-24-13.0 

2.3283566 

-3.4999056 

-0.50715805 

2.2721038 

1.2746548 

1.6158241 

Arc 3( 

64-4-1-22.0 

0.76296064 

3,6490419 

1.9857054 

-2.7806607 

1,0205339 

-0.81923611 


64-4-7-15.0 

3.4461234 

-2.4088100 

0.50885287 

1.3277141 

2.2415368 

1,6144312 


1^64-4-13-19.0 

3.38 7 2421 

1.3385265 

2.1522760 

-1.4069923 

2.6725021 

0.54080057 


44-4-25-2.0 

-2.5603443 

3.3390780 

0.29994066 

-2.1318214 

-1.4887792 

-1.6428786 


64-4-28-15.0 

3.2628949 

-2.6693162 

0,32101150 

1.5457882 

2.0856009 

1.6366277 


64-5-5-16.0 

3.7257132 

-1.7632128 

0.94048572 

0.82412264 

2.5449717 

1.5038900 


64-5-12-16.0 

3.7475309 

-1.6867283 

0.9 9 296812 

0.76273664 

2.5725322 

1.4891387 

Arc 4( 

64-5-19-14.0 

2.7849331 

-3.1784434 

-0.10945226 

1.9747957 

1.6719452 

1 .6505955 

1 

64-5-25-15.0 

3.4307022 

-2.4 1 26 224 

0.53392194 

1.3319460 

2.2493385 

1.6072617 


64-6-2-21.0 

1.7406646 

3.1597779 

2.1946670 

-2.5288299 

1.6899106 

-0.43439828 


64-6-9-21.0 

1.6500071 

3.2194181 

2.1774956 ■ 

-2.5601432 

1.6304633 

-0.47782909 


64-6-16-15.0 

3.5655116 

-2.1536100 

0.72108034 

1.1248629 

2.3871212 

1.5668312 


1^4-6-23-15.0 

3.6075927 

-2.0587512 

0.78733380 

1.0494437 

2.4320607 

1.54<15642 


'44-7-4-2.0 

-3.2093967 

2.7165420 

-0.32063943 

-1.5841089 

-2.066085U 

-1,6378223 


64-7-7-3.0 

-3.6954496 

1.7950587 

-0.94983861 

-0.85046154 

-2,5458831 

-1.5023067 


64-7-13-17.0 

3.7354217 

0.41539202 

1.9259618 

-0.76727529 

2.8431082 

0.87209834 


64-7-21-21.0 

0.68011374 

3.6978808 

1.9157389 

-2.7798310 

-0,96706741 

-0.88586596 


64-7-27-16.0 

3.8417868 

-0.11824881 

1.7510381 

-0.39126192 

2.8593634 

1.0493355 


64-8-3-17.0 

3.4664817 

1.1480683 

2.1210692 

-1.2649191 

2.7375684 

0,58214636 


64-8-11-1.0 

-3.3701691 

2.4805650 

-0.51870386 

-1.3894809 

-2.2250586 

-1.6062097 


64-8-17-19.0 

1.6959231 

3.2024028 

2.1645995 

-2.5377008 

1.6647695 

-0.48003497 


64-8-25-10.0 

1,5330007 

-3.8258893 

-0.91436182 

2.5905352 

0.67495556 

1.5073075 


64-9-1-10.0 

1.7622665 

-3.7590544 

-0.76800569 

2.5108491 

0.85857627 

1.5478150 


64-9-9-14.0 

3.8729406 

-0.56806750 

1.5864037 

-0.064200095 

2.8378419 

1.1715632 


64-9-15-12.0 

3.4973065 

-2.2642849 

0.69291203 

1.2114165 

2.3483095 

1.5648130 


64-9-22-10.0 

2.3957332 

-3.4619283 

-0.32713684 

2.2142000 

1.3764322 

1.6238887 

Arc 5/ 

64-9-29-6.0 

-1.2057564 

-3.4918173 

-2.0387967 

2.6768512 

-1.3377786 

0.70292893 


64-10-6-5.0 

-1.9027758 

-3.0642587 

-2.1871975 

2,4522966 

-1.8112966 

0.39976832 


64-10-13.0 

-3.7909750 

1,3890374 

-1.2159474 

-0.54220619 

-2.6929405 

-1.3843106 


64-10-20-16.0 

2.3418688 

2.7061755 

2.2406431 

-2.2408225 

2.0926322 

-0.18928830 


64-10-26-16.0 

2.1737218 

2.8588045 

2,2200917 

-2.3278640 

1.9849102 

-0.28044993 


64-11-2-5.0 

-1.0575392 

-3.5703165 

-1.9843669 

2.7067256 

-1.2347065 

0.77366371 


64-11-11-2.0 

-3.2269498 

-1.5931919 

-2.1991397 

1.5433675 

-2.6360297 

-0.35716605 


64-11-17-6.0 

0.44169064 

-3.9493268 

-1.4227193 

2.7947443 

-0.14738721 

1.2691720 


65-1-10-6.0 

1.8733482 

-3.7329199 

-0.61394837 

2.4551255 

0.97133874 

1.5716495 


65-1-13-16.0 

-0.024904471 

3.9169200 

1.5734233 

-2.8064159 

0.44943751 

-1.1686893 


65-1-20-12.0 

3.2361444 

1.5906403 

2.1989944 

-1.5260383 

2.6448765 

0.33028481 


65-1-27-4-5.0 

0.42866892 

-3.9652039 

-1.3910958 

2,7914381 

-0.14522612 

1.2721610 


65-2-2-13-30.0** 

I 1.8326846 

3.1377111 

2.1684747 

-2.4504221 

1,7730384 

-0.49949115 


v6 5-2-16-4-5.0 

j 0,94696655 

-3.9573969 

-1.1249388 

2.7174166 

0.24662704 

1.4128023 


'From thr Data Systems and Tracking Directorate, computed from range and range rate and Minitrack Data by Robert Chaplick, Gerald Repass, and Catleton Carver, Iron an orbit determination program 
due to Dr. Joseph Siry (see text in Appendix A of Reference 14 j for earth- gravity constants used in this program.) The iiieccial orthogonal system x, y, z has the z axis pointing towards the vernal 
equinox of epoch and the z axis pointing towards the North Celestial Pole. 

’*This orbit gave unacceptably large residuals in the acceleration analysis in Section 1 and was consequently ignored there in the final reduction of arc ^ data. 
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Table A2 

Inertial Position and Velocity Coordinates or Mean Elements for Syncom 2 and 

Syncom 3 as Reported by GSFC* 
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Table A3 

Mean Elements and Related Data for Syncom 2 in Arc 8, As Reported By GSFC 


1 (Descending Eq. Crossing Data) 

Time From 
65-5-3-13-27- 
15.5:Sid. Days 

-28 

-21 

-14 

-7 

0 

7 

14 

21 

28 

35 

Time; 

(Yr-Mo-Day- 
Hr- Min- Sec, U.T.) 

65-4-5-15-16-37.1 

65-4-12-14-50-12.8 

65-4-19-14-21-50.3 

65-4-26-13-54-55.1 

65-5-3-13-27-15.5 

65-5-10-12-59-47.7 

65-5-17-12-32-26.3 

65-5-24-12-4-40.5 

65-5-31-11-8-13.5 

65-6-7-11-8-13.5 

[Descend- 
ing Eq. 
Cross.) 
Deg's. 

66.023 

65.794 

65.647 

65.4895 

65.362 

65-265 

65.200 

65.141 

65.1575 

65.183 

Eccentri- 

city 

e 

(10-5) 

63 

53 

51 

56 

R.A. of The 
Asc. Node 
0 

(Deg’s) 

308.312 

308.192 

307.954 

307.882 

Arg. of 
Perigee 

<o 

(Deg’s) 

336.214 

345.151 

334.768 

337.312 

Mean 

Anomaly 

MA 

(Deg’s) 

15.179 

13.186 

30.730 

35.182 

Inclination 

i 

(Deg’s) 

31.956 

.930 

.939 

.912 

.955 

.830 

.848 

.867 

.787 

.807 

.731 

31.820 

.717 

.725 

.727 

Avg: Orbits 
8-6/8-15; 
i = 31.786 
Avg: Orbits 
8-7/8-15 
i = 31.780 

(Asc. Eq. Crossing Data) 

X 

(Ascending 
Eq. Cross.) 
Deg’s 

67.733 

.423 

.270 

.007 

66.314 

.071 

65.822 

.712 

.534 

.404 

.288 

.244 

.155 

.1895 

.217 

Time From 
Jan. 109.0997, 
1965 (Sid. Days) 

-53 

-47 

-44 

-37 

-21 

-14 

-7 

0 

7 

14 

21 

28 

35 

43 

50 

Time: 

From 

1965.0 

(Days) 

56.24 

56.2412 

62.2257 

65.2177 

72.1986 

88.1567 

95.1378 

102.1195 

109.0997 

116.0811 

123.0619 

130.0428 

137.0238 

144.0046 

151.9820 

158.9626 

Estimated 
(10"® E.R.) 

1.650 

.755 

13.3 

1.680 

7.22 

1.45 

1.70 

0.567 

.810 

4.51 

1.63 

2.97 

1.67 

7.66 

2.09 

Semimajor 

Axis 

(a): 

Earth Radii 

6.6110989 

3906 

6.6109100 

6.6112847 

3429 

1906 

1446 

1335 

0585 

0267 

0035 

6.6109407 

8998 

8153 

7771 

7276 

Avg: Orbits 
8-6/8-15: 
a = 6.610953 
Avg: Orbits 
8-7/8-15 
a = 6.610931 

Tracking 

Epoch; 

(Yr-Mo-Day- 
Hr-Min: U.T.) 

65-2-24-19-30.0 

65-2-25.0 

65-3-3.0 

65-3-6.0 

65-3-13.0 

65-3-29.0 

65-4-5.0 

65-4-12.0 

65-4-19.0 

65-4-26.0 

65-5-3.0 

65-5-10.0 

65-5-17.0 

65-5-24.0 

65-5-31.0 

65-6-7.0 

Orbit, 

8- 

1-H I-H CQ CO m CD D»00 O>Oi-HCQCO^IA 


§ 





•o 



3 

2 
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Arc 8-2 Analyses the Descending Equator Crossing data from orbits 8-6 through 8-15 (slow drift) 




Table A4 

Mean Elements And Related Data for Syncoms 2 and 3 In Arcs 1, 

as Reported By GSFC 


Mean a 
(E.R.) 

6.6107434 

9495 

462 

6.6111018 

0200 

763 

6.6113726 
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130 

P7Q 

6.6118417 
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0819 

3652 

3870 
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Nov. 1963* 
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00 csi o w oj «D o eo IT) od o CX5 eo 

^^eococ^i 1 1 

1 1 1 1 1 

Epoch 
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Orbit, 

2- 

i-tc<lco^in«?c-ooo>Oi-HC>3CO'^msc> 


Time 
From Jan. 
329.6393,1964 
(Sid. Days) 

'^i-Ht-ooOcoirtC'lt- 
CM CM 1 tH CM CM 

1 1 i 

Long, of 
First 
Asc. Eq. 
Cross 
After 
Epoch: 
k 

(Deg's) 

180.219 

.028 

179.737 

.191 

178.798 

.509 

.171 

177.968 

.814 

Time* 

From 1964.0, 
(Days) 

305.9723 

8.7567 

12.8308 

21.6437 

29.6393 

35.5028 

44.4400 

51.4760 

356.4394 
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mated 

cr(a): 
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E.R.) 
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CM 
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a 

(E.a) 
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Orbit, 

6- 
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From 1965.0 
(Days) 
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47.7895 

54.8054 

61.732 

68.8539 

75.7421 

ta o 

w 

2.410 

1.150 

1.330 

0.453 

0.561 

0.682 

0.473 

0.495 

Mean 

a 

(E.R.) 
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Long, of 
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Asc. Eq. 
Cross 
Past 
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(Deg's) 
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-158.96 
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Mean a 
(E.R.) 

6.6105202/6.6104741 

05391 

298 

7646 

979 

8160 

472 

702 

9835 

6.6112494/6.6109979 

126 

543 

3540 

6250 

5592 

7250 

815 

6.6120076 

Time in 
E>ays From 
63-8-68.495* 

-50.365 

-45.379 

-41.390 

-37.401 

-32.415 

-28.427 

-24.438 

-19.451 

-16.459 

-9.478 

-5.489 
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7.477 

15.455 

23.435 

30.416 
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42.386 

Epoch 

(Yr-Mo-Day- 
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♦Time at first ascending equator crossing past the epoch. 






Table A5 


Mean Elements and Related Data for Syncom 2, Arc 5, 
as Reported by GSFC 


Orbit # 
5- 

Tracking 
Epoch 
(See Ref. 
[4]) 

Mean 

Semimajor 
Axis, a : 
(Earth Radii ^ 

. t 

\ (VSid.Day) 

i 

Inclination 

(Rad.) 

\ 

(First Asc. 
Eq. Cross. 
Past Epoch) 
Deg's 

Est. 

cr(a): 

(10"5 E.R.) 

. t 

(ySol.Day) 

1 

64-7-4-2.0 

6.6165547 

-0.48480 

.56794437 

-171.26 

.5440 

-.4861 

2 

64-7-7-3.0 

6.6166962 

9636 

.56787141 

-172.70 

.0920 

977 

3 

64-7-13-17.0 

6.6165633 

8551 

.56726514 

-176.04 

.0418 

868 

4 

64-7-21-21.0 

4925 

7972 

50460 

-179.79 

.1040 

810 

5 

64-7-27-16.0 

4312 

7471 

26786 

177.41 

.0610 

760 

6 

64-8-3-17.0 

3966 

7189 

.56681234- 

174.18 

.0623 

732 

7 

64-8-11-1.0 

3059 

6448 

.56628775 

171.00 

.0683 

658 

8 

64-8-17-19.0 

3334 

6673 

563818 

167.86 

.105 

680 

9 

64-8-25-10.0 
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164.69 

.370 

631 

10 

64-9-1-10.0 
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543507 

161.57 

.129 

665 

11 
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157.51 

.168 

558 
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6632 
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154.82 

.118 

676 
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152.09 

.179 
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148.87 

.140 
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145.70 

.152 
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72.25 

1.260 

125 


Avg: 6.616855 


Avg: 32.33“ 





Use 

: a = 


Use i= 32.3“ 






6.6169 







^ Estimated from Mean a data. 
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Table A 6 



172.558 



Table A 7 

Syncom Vector Elements as Reported by DOD, Air Force Systems Command 
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*M = Month, D = Day Y = Year, H = Hour. MI = Minute. S = Second. 
Also, see notes at Bottom ol Table A8 




Table A8 

Syncom Vector Elements as Reported by DOD, Air Force Systems Command 


Syncom 3 

1 Date (U.T.) 

a" 



AZ" 

RAD' 

V/FPS 

M 

D 

Y 

H 

MI 

S 

2 

4 

65 

0 

0 

0 

J329.3087100 

.672138 

89.9941822 

89.923047 

138375573.3 

10085.61611 

3 

5 

65 

0 

0 

0 

337.2338772 

.0907950 

89.9957186 

89.9519873 

138387194.8 

10085.00939 

3 

11 

65 

0 

0 

0 

342.465331 

.0930934 

89.9927339 

89.9545506 

138387868.4 

10084.73452 

3 

6 

65 

0 

0 

0 

346.9072271 

.3243960 

90.0036182 

90.040089 

138442281.5 

10081.23290 

3 

24 

65 

0 

0 

0 

353.5587907 

.0634430 

89.9974880 

90.0087157 

138364442.3 

10085.98775 

3 

31 

65 

0 

0 

0 

-.0297323 

-.0466264 

89.9945489 

90.0218809 

138309091.8 

10089.18472 

3 

31 

65 

0 

0 

0 

225.6557090 

-.0349054 

90.0046737 

89.9631400 

138318082.0 

10088.34542 

7 

8 

65 

5 

11 

0 

173.8584686 

.1177573 

90.0007217 

89.9671938 

138333100.1 

10087.65142 

7 

14 

65 

13 

25 

0 

303.6761356 

-.0573860 

89.9978681 

90.1236198 

138335567.3 

10087.64151 

7 

20 

65 

13 

25 

0 

309.7154355 

-.0620873 

89.9971435 

90.1190876 

138336828.1 

10087.21652 

9 

2 

65 

14 

25 

0.0 

9.7560703 

-.2697071 

90.0000097 

90.0123450 

138335712.4 

10087.11885 

9 

29 

65 

13 

41 

0.0 

26.9170975 

-.2961765 

90.0012623 

89.9495137 

138334916.4 

10087.08653 

10 

15 

65 

1 

30 

8.883 

219.6767686 

.2590034 

89.9986513 

90.0395497 

138280569.8 

10091.82446 

11 

11 

65 

0 

1 

14.875 

222.2051107 

.2747621 

90.0001025 

90.0741580 

138312827.4 

10090.10942 

11 

23 

65 

23 

32 

17.739 

226.8336731 

.2761059 

90.0003455 

90.1791327 

138316354.7 

10089.77591 

12 

9 

65 

23 

25 

6.427 

239.7850209 

.2948162 

90.0039197 

90.3207487 

132315950.9 

10089.68714 

Corrected Orbit 


(6.384) 

(239.7821569) 

(.2834610) 

(90.0044185) 

(90.3277418) 

(138315042.5) 

(10089.75657) 

12 

21 

65 

23 

0 

8.657 

244.6624456 

.2502542 

90.0005093 

90.3725080 

138317347.2 

10059.55986 


(22) 

Cor. 

Orb. 

(11) 

(46.219) 

(248.5289543) 

(.2389576) 

(90.0009380) 

(90.3940341) 

(138319005.7) 

(10089.41514) 

1 

5 

66 

22 

59 

^5.906 

258.5706397 

-.0085589 

89.9936737 

90.2325447 

138287344.2 

10091.68830 

C. Orb 

.(6) 


(23) 

(4) 

(31.968) 

(260.7895521) 

(-.0158445) 

(89.9935676) 

(90.2373577) 

(138288609,2) 

(10091.57113) 

1 

19 

66 

20 

56 

16.662 

240.8567046 

.0717809 

90.0022661 

90.2796546 

138295615.7 

10090.89758 

2 

4 

66 

20 

44 

20.747 

252.9875511 

.2019325 

90.0069535 

90.5668549 

138312648.7 

10089.54462 

3 

3 

66 

22 

1 

14.254 

298.0440906 

-.2997467 

90.0026493 

90.6568832 

138310811.7 

10089.77604 

3 

17 

66 

22 

34 

35.662 

319.8424334 

-.5453541 

90.0018012. 

90.56183511 

138319644.1 

10089.08818 

3 

30 

66 

21 

24 

23.764 

314.7805399 

-.5317386 

89.9994567 

90.5276043 

138313388.1 

10089.48471 

4 

10 

66 

0 

13 

2.374 

6.7218328 

-.7224497 

90.0011987 

89.9537090 

138329513.3 

10088.16882 

5 

2 

66 

23 

28 

43.618 

17.9376586 

-.7580983 

90.0018212 

89.8618606 

138324944.6 

10088.64166 

5 

17 

66 

0 

1 

44.652 

39.8916786 

-.7961832 

89.9983801 

89.4885543 

138318036.3 

10089.08628 

6 

1 

66 

0 

56 

32.305 

68.2949586 

-.4778375 

89.9980245 

89,1712197 

138310369.7 

10089.60206 

6 

20 

66 

1 

8 

18.168 

89.6999357 

-.0837705 

89.9974096 

89.0012170 

138320870.1 

10089.34540 

6 

30 

66 

23 

51 

25.402 

80.6468849 

-.2472802 

90.0014671 

89,0450739 

138312020.7 

10089.94653 

7 

14 

66 

17 

37 

40.716 

0 

-.9924101 

90.0129977 

90.1462207 

138353051.0 

10086.69946 

9 

16 

66 

18 

5 

32.019 

67.8130270 

-.6941001 

90.0027906 

88.9776683 

138317331.3 

10083.92757 

10 

14 

66 

11 

45 

58:570 

0.0 

-1.2705583 

90.0089077 

90.1510399 

138337220.9 

10087.68599 

12 

13 

66 

IS 

40 

21.168 

117.3088697 

.5734308 

89.9964678 

88.6843516 

138334769.4 

10087.66824 

1 

1 

67 

0 

1 

0.0 

260.3029154 

.3491017 

90.0037591 

91.4563330 

138340290.5 

10087.30701 


♦The elements refer to the equinox and equator of date. 

The earth model used to determine the elements, from data over 
about a week's duration (maximum), has the following parameters: 


fx = .0053045314 E.R.Vmin^ 

Ro = 20925738.19 FT/E.R. 

I Km.= 3280.8399 Ft. 

Zonal gravity field through Jg: Ref., W.H. Guier (1965), "Recent Progress in Satellite Geodesy" Applied Physics Lhb. Report 
TG-659, Feb. 1965, Johns Hopkins Univ., Silver Spring, Maryland. 

Non zonal gravity field through J44: Ref. (Prior to spring 1966); W.H. Guier & R.R. Newton (1963), "Non Zonal Harmonic 
Coefficients of the Geopotential from Satellite Doppler Data", APL John's Hopkins Univ. Report # TG-520 (1963) Ref., (after 
spring 1966); Guier & Newton (1963) except for H22 and H33 Harmonics which are from C.A. Wagner (1966), "Longitude 
Variations of the Earth's Gravity Field as Sensed by the Drift of Three Synchronous Satellites", Journal of Geophys. Res. , 
71, #6, 1703 (1966). 

Non zonal gravity field from J51 to JgB: Ref., W.H. Guier (1965), used after spring 1966. 
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Syncom 2 Orbit Data in Arc 8 as Reported by DOD 
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♦Right ascension of the ascending node: **Hour angle of Greenwich: ( ) = Estimated, not reported, data 
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Table All 

Syncom 3 Orbit Data in Arc 13, as Reported by DOD 


Date 

Longitude of 
Asc. Eq. 
Cross 
(Deg's) 

Time 
From 6 
Sept. 1966 
(Days) 

Observed 
From Eq. 
Cross. Long. 
Shift in One 
Revolution 
(Deg./Day) 


8 Sept. 1966 

161.28 

2 

-0.005 



(161.245) 

(5.5) 


^-0.010) 

15 Sept, 1966 

161.21 

9 

-0.010 



(161.15) 

(16) 


(-0.0085) 

29 Sept. 1966 

161.09 

23 

-0.009 



(161.02) 

(30.5) 


(-0.0095) 

14 Oct. 1966 

160.95 

38 

-0.013 



(160.655) 

(60) 


(-0.010) 

13 Dec. 1966 

160.36 

98 

-0.012 



(160.275) 

(105) 


(-0.012) 

27 Dec. 1966 

160.19 

112 

-0.012 



(160.065) 

(119.5) 


(-0.0165) 


11 Jan. 1967 


159.94 


127 


-0.016 



Table A12 

Equator Crossing Data and Sets of Elements for Early Bird in Arcs 9 and 12 
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Table A12 (Cent.) 

Equator Crossing Data and Sets of Elements for Early Bird in Arcs 9 and 12 
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Table A12 (Cent.) 

Equator Crossing Data and Sets of Elements for Early Bird in Arcs 9 and 12 
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APPENDIX B 

Estimating the Harmonic Variability 


From the covariance statistics, derived as a by-product of the least squares 
fitting process, we can estimate the degree of harmonic variability which the data 
allows. In this calculation, the model is assmned to be unbiased with respect to 
the data. The data residuals are assumed to be uncorrelated and random normally 
distributed. In addition, the measurement variances are assumed to produce 
realistic weights for the data in those model tests where weighting is employed. 
Actually, we know the abbreviated models tested cannot yield unbiased residuals 
with this limited data, or uncorrelated residuals, due to the effects of the har- 
monics ignored in the fit. Nevertheless, we can hope to learn something of the 
variability in the abbreviated model, allowed by the data, if we proceed in two 
stages. In the first stage, we assume that all the unmodeled effects are zero. 

If this were the case, our assumptions as to unbiasness, noncorrelation and 
realistic weights would be true. In the second stage of assessing the modeled 
harmonics variability, we can actually account for the unmodeled effects by 
using likely values for them from outside sources. If these outside-source values 
are reasonably realistic, then again the residuals should be unbiased and uncor- 
related, and the statistics of the solution would be valid. 

In this report, the covariance statistics are chosen from the solution without 
any assumption as to unmodeled effects (Solution 6, Table 2; also see Figure 2). 
(The general least squares program used, and the derivation of the statistical 
data from it, is described in Reference [Bl] .) In this way, we maintain maximum 
faithfulness to the actual data; for example, in the determination of accelerations of 
new 24 hour satellites, and in comparison of the results with other determinations 
from different data. However, we acknowledge the likely contribution of unmodeled 
effects by observing the gross changes possible in the modeled harmonics, when 
reasonable values for these effects are included in least squares solutions. This 
second stage process has been described already in the text. The best harmonic 
solution (at the bottom of Table 2) has bounds which are set by considering both 
the internal variability of the solution (first stage error analysis) and the varia- 
bility of the harmonics themselves in tests with reasonable account taken of un- 
modeled, higher order resonance effects. 

A more desirable and logical way to conduct the first stage error analysis 
(no specific assumptions on higher order terms) would be to increase all the 
data deviations by an expected RMS (root mean square) deviation due to the 
unmodeled effects. While this would not eliminate the bias and correlation of 
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the residuals in the first stage analyses, it would result in more realistic weights 
to the data and a better overall solution. For example, Reference [4 ] shows that 
the unmodeled 4th order harmonics may produce accelerations as high as .016 x 
10“® radians/day2 (H^^) and .018 x io~® rad./day^ (H^j) on the geostationary 
satellite. In gravity solutions only for and , an expected data devia- 

tion (RMS), due to and , of .707 (.OIG^ + .OIS^)!/^ x lO"® = .017 x lO'^ 
rad/day ^ should be allowed for in the measurements on the geostationary 
satellites. Measurements on the moderately inclined Syncom 2 in this limited 
model test need allow for only .009 x 10“® rad/day^ RMS deviation due to 
and principally because the effect of i® almost negligible at the inclina- 
tion of that satellite. Future 24 hour satellite limited gravity tests will be made 
with increased data variance to achieve better overall solutions. 


1(7 Ellipses 

The li7 ellipses in Figure 2 summarize the covariance statistics of the 
basic (six parameter, three gravity term) solution (first stage error analysis) as 
it applies to the permitted variability of each harmonic term as seen by the data. 
Each harmonic term is specified by the coefficients and , so that the 
statistics necessary to define its variability are: 

Cd , Sd , crC j> , (7Sp andc(Cp , Sp ); 

'Vm 'Vm 'om 

the first two symbols being the least squares estimates of the harmonic coef- 
ficients, the second two, their standard deviations and the last, the correlation 
coefficient (-1 < c < 1). Providing the experiment residuals are random 
normally distributed, it is shown in Reference [B2 ] that the above statistics 
define a probability ellipse, X , in the C p , Sp plane defined by: 

'wxu 


+ 




_^2 


(Bl) 
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This ellipse has the property that on it, the probability density function for 
possible sets of and seen by the data, is constant. The ellipse scale is 
set by which defines the probability P (7) of finding sets of and with- 
in the ellipse. It is shown in Reference [B2 ] by a simple integration of the 
density function that this probability P(7) is related to 7^ and c by the formula: 


P(7) = 1 - exp [- 7 V 2 (l-c2)] . 


(B2) 


Thus, for example, the la probability ellipse contains an area in which it is ex- 
pected that 68.26% of the sets of and values will fall. Evaluating (B2) 
for P (7) = 0.6826, 

72 ^ = 2.30 (l-c2). (B3) 


Similarly, the 2a ellipse [P(7) = .9544] defines the scale size from: 

72 ^ = 6.16 (l-c2) (B4) 

A rather simple exercise in analytic geometry shows that the probability 
ellipse of Equation (Bl), with respect to the plane: 

1) Is centered at Co , So , 

't/m 'Vm 

2) Has major (A) and Minor (B) axes given by: 


A = 27V {D- [E2 + 4c2/f2]V/2} 


(B5) 


where 


B = 27V{D+ [£2 + 4c2/p2] } 


1 

1 

+ . . 



1 

1 




(B6) 


(B7) 
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and 


F = 


and 

3) Has it*s minor axis (B) rotated counterclockwise from the axis by an 
angle 6 , given by: 


6 


- TAN-i 
2 


-2c/F 

E 


(B8) 


These formulas have been used in constructing the probability ellipses in Figure 2 
from the covariance statistics of the basic solution (6) in Table 2. 

After considering the likely model bias effects (in the second stage of the 
error analysis) on the basic solution, somewhat subjective upper and lower 
realistic bounds have been set on the modeled coefficients. They are displayed 
by the dashed boxes in Figure 2. Defining the absolute variability possible in 
by the boxed bounds, does not take into account the correlation between 
and in the various tests used to set these bounds. Therefore, the 
boxed bounds for should be reasonably conservative since the set of , 

S^in values used to define it do not (in general) tend to fill the inner space in a 
random manner. To insure this conservatism in reporting the and 
coefficients (which are the polar coordinates of points in the and plane), 

I have calculated the possible extremes of their variation by considering rays 
to the corners of the boxed bounds in Figure 2. In future reports, the correlation 
of coefficient solutions will be taken into account when the unmodeled effects are 
estimated. When this is done, the absolute bounds on the gravity solution should 
be known more precisely. 
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APPENDIX C 
List of Symbols 

Specifying the amplitude and phase angle of the spherical 

gravitational harmonic whose order (or latitude frequency) is 

■I and whose longitude frequency is m. 

Referring to the spherical gravitational harmonic of order ^ and 
longitude frequency m . 


dx d^x 
dt ’ 

Semimajor axis of the satellite's orbit, synchronous semimajor 
axis (a^ = 6.611 earth radii). 

Inclination of the satellite's orbit. 

Standard deviation. 

Velocity increment. 

Parameter defining the scale of the probability ellipse. 
Correlation coefficient. 

Earth radii. 

Root mean square. 

Ascending equator crossing. 

Descending equator crossing. 
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